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The Photon Dominated Region in the IC348 molecular cloud 
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ABSTRACT 

Aims. In this paper we discuss the physical conditions of clumpy nature in the IC 348 molecular cloud. We show that the millimeter and 
sub-millimeter line emission from the IC 348 molecular cloud can be modelled as originating from a photon dominated region (PDR). 
Methods. We combine new observations of fully sampled maps in [C i] at 492 GHz and 12 CO 4-3, taken with the KOSMA 3 m telescope 
at about 1' resolution, with FCRAO data of I2 CO 1-0, l3 CO 1-0 and far-infrared continuum data observed by HIRES/IRAS. To derive the 
physical parameters of the region we analyze the line ratios of [Ci] 3 P 1 - 3 P / 12 CO 4-3, [Ci] 3 P,- 3 P / 13 CO 1-0, and 12 CO 4-3/ 12 CO 1-0. A 
first rough estimate of abundance is obtained from an LTE analysis. To understand the [C i] and CO emission from the PDRs in IC 348, we 
use a clumpy PDR model. With an ensemble of identical clumps, we constrain the total mass from the observed absolute intensities. Then 
we apply a more realistic clump distribution model with a power law index of 1.8 for clump-mass spectrum and a power law index of 2.3 for 
mass-size relation. 

Results. We provide detailed fits to observations at seven representative positions in the cloud, revealing clump densities between 4 10 4 cirT 3 
and 4 10 s cirr 3 and C/CO column density ratios between 0.02 and 0.26. The derived FUV flux from the model fit is consistent with the field 
calculated from FIR continuum data, varying between 2 and 100 Draine units across the cloud. We find that both an ensemble of identical 
clumps and an ensemble with a power law clump mass distribution produce line intensities which are in good agreement (within a factor ~ 2) 
with the observed intensities. The models confirm the anti-correlation between the C/CO abundance ratio and the hydrogen column density 
found in many regions. 
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1. Introduction 

Photon Dominated Regions (PDRs) are surfaces of molec- 
ular clouds where chemistry and heating are regulated by 
far-ultraviolet (FUV) p hoton s (6.0 eV < hv < 13.6 eV ) 



dHollenbach & Tielensl 1 1999b . Gas cooling happens mainly 



via the fine struct ure lines of [Oil, [Cn l, [Ci] and the rota- 
tional lines of CO (Kaufman et all 1 19991) . The form of carbon 
changes with increasing depth from the surface of the PDR 
from C + through C° to CO. Therefore emission from [C n], [C i] 
and the rotational lines of CO can be used as pro bes of temper- 



ature, density and column de nsity in the PDRs (IRramer et al 



2004t lMookeriea et all [20061) . In particular, the two fine struc- 



ture lines of neutral atomic carbon at 492 and 809 GHz give in- 
formation on intermediate regions between atomic and molec- 
ular gas, since they trace both physical regimes. Initial ho- 
mogeneous plane-parallel steady-state models of PDRs pre- 
dicted that PDRs are confined to the FUV irradiated surfaces of 
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molecular clouds. However, observations detected widespread 
[C n] and [C i] emission, thus suggesting that PDR surfaces are 
formed deeply inside the molecular clouds. One possible ex- 
planation is in terms of the dumpiness of the molecular clouds 
which allows for deeper penetration of FUV radiation so that 
PDRs are formed at the surfaces of the clumps present through 



198alHoweetallll99U 



out the molecular cloud dStutzki et al 

Mookerjea et all 120061) . In addition, an explanation has also 

been proposed in terms of the dynamic s of the ISM: a. clo uds 



are suddenly shielded from FUV fields dStorzer et al.Lll997l) : b. 



time-dependent chemistry of the PDRs (IQka et all 120041) . 
IC348 is a young open cluster northeast of the 



Perseus cloud at a distance of 316pc (iHerbigl 119981) . The 
IC 348 cluster has about 300 stellar members, with spec- 
tral types fro m B5 to late M, as found from observa- 
tions in X-rays (IPreibisch et all 1996t Preibisch & Zinnecker 
200lh. near- i nfrared dLada & Lada[ll995l:lLuhman et a J J1998 



Naiita et all l2000l: JLuhman et all l2005l). a nd op tical (IHerbig . 
119981: iLuhmanetall 120031) . iLadaeTall d2006l) constructed 
optical-infrared spectral energy distributions (SEDs) for all 
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known members of the cluster by combining the Spitzer Space 
Telescope and ground-based observations. The IC348 molec- 
ular cloud prominently shows up i n the CO 3-2 su rvey of the 
Perseus molecular cloud complex (ISun et all 120061) . The FUV 
radiation of its youngest member, HD 28 1159, a B5 V star, 
heats the surrounding gas and dust. Surveys for molecular out- 
flows revealed so far four bipolar outflow s powered by deeply 



embe d ded young sources (|Eisloffel et all l2003t iTafalla et al„ 
2OO6l) . |j0rgensen et alJd2007l) estimated the current star forma 



tion efficiency of 10%— 15% for the whole Perseus cloud using 
Spitzer mid-infrared data and SCUBA sub-millimeter maps. 



Sun et alj ( 120061) observed IC 348 in four low-/ CO tran- 



sitions as part of their 7.1 square degree survey of the entire 
Perseus molecular cloud and quantified the spatial structure of 
the emission in t erms of its power sp ectrum, applying the A- 
variance method dStutzki et all Il998|) to the integrated inten- 
sity maps. The A-variance spectra of the overall region was 
found to follow a power law with an index of /3 - 2.9 - 3.0; 
while sub-regions of Perseus gave different power spectral in- 
dices. The active star forming regions showed higher indices 
typical for large condensations; lower indices were derived for 
dark clouds, representing more filamentary structure. IC 348 as 
an exception, showed a rather filamentary structure correspond- 
ing to the smallest value of/? = 2.71 in 13 CO 2-1. 

In this paper we present large-scale (20'x 20') fully- 
sampled maps of the [Ci] 3 Pi - 3 P (hereafter [Ci]) and 12 CO 
4-3 emission from the IC 348 molecular cloud. At a common 
resolution of 70", we have combined our data with the 12 CO 
1-0, 13 CO 1-0 data fr om the Five Colleg e Radio Astronomy 
Observatory (FCRAO) JRidge et alll2003l) . and the far-infrared 
(FIR) continuum data from HIRES/IRAS. The goal of the anal- 
ysis is to understand the extent to which the observed line emis- 
sion from IC 348 can be understood in the framework of a pho- 
ton dominated region. 

Section|2]describes the KOSMA observations and the com- 
plementary datasets used for analysis. Section [3] discusses 
the general observational results. A Local Thermodynamic 
Equilibrium (LTE) analysis of the observed line intensities is 
presented in Section |U A comparison with clumpy PDR mod- 
els is presented in Section[5] Section[6]summarizes the results. 

2. Datasets 

2.1. [C \] and n CO 4-3 observations with KOSMA 

We have used the Kolner Observatorium filr Sub-Millimeter 
Astronomie (KOSMA) 3-m sub-millimeter t el escope 
on Gornerg rat, 
Kramer et al. . 



Switzerland 



(IWinnewisser et al. , 
1998a) to observe the emission of the fine 



mm 



structure line of neutral carbon at 492 GHz and 12 CO 4-3 at 
461 GHz. We have used the Sub-Millimeter Array Receiver for 
Two frequencies ( SMART) on KOSMA for these observations 



dGraf et all |2002|) . SMART is an eight-pixel dual-frequency 



SIS-heterodyne receiver capable of observing in th e 650 



and 350 yum atmospheric windows dGraf et all 12002b . The 



IF signals were analyzed with two 4x 1 GHz array-acousto- 
optical spectrome ters with a spectral resolution of 1.5 MHz 
dHorn et all 1 19991) . The typical double side band receiver noise 



temperature at 492 GHz is about 150 K. The observations were 
performed in position-switched On-The-Fly (OTF) between 
December 2004 and February 2005. Owing to technical 
difficulties the higher frequency channel of SMART could not 
be used at the time of these observations. 

We observed a fully sampled map of the IC 348 molecular 
cloud centered at a = 03 /, 44 m 10' 5 , 6 = 32°06'(J2000), extend- 
ing over 20' x 20'. For the observ ations, we used the p osition- 
switched on-the-fly (OTF) mode dBeuther et all 120001) with a 
sampling of 29". The emission-free off position was selected 
from the 13 CO 1-0 FCRAO map, which is (-8',10') relative to 
the map center. We estimate the pointing accuracy to be bet- 
ter than 10". The half power beam width (HPBW) and the 
main beam efficiency (B e ff) were derived from continuum scans 
of the Sun and the Jupiter. The HPBW at both frequencies is 
60"and the main beam efficiency is 50%. The forward effi- 
ciency F e ff of the telescope is ~ 90%. Atmospheric calibration 
was done by measuring the atmospheric emission at the OFF- 
position and using a standard atmospheric model to fit the opac- 
ity taking into accou nt the sideband imbalance dCernicharol 
1985HHivam"aill998b . 

All data presented in this paper are in units of main beam 
temperature r m b, calculated from the observed calibrated an- 
tenna temperature T^ using the derived beam and forward effi- 
ciencies, r m b = T* k (Feff/Seff)- Based on observations of refer- 
ence sources such as W 3 and DR 21 we estimate the accuracy 
of the absolute intensity calibration to be better than 15%. The 
data reduction was carried out using the GILD AS Lj software 
package. 

2.2. Complementary data sets 

We have used the FCRAO 12 CO 1-0 and 13 CO 1-0 datasets 



with an resolution of 46"presented by iRidge et al 



comparison. For comparison, we have smoothed all CO and 
C i data to a common resolution of 70"assuming that both the 
FCRAO and KOSMA telescopes have a Gaussian beam. 

Further, we have obtained HIRES processe d IRAS maps of 
the du st continuum emission at 60 and 100/zm dAumann et all 
19901) . The dust continuum maps have resolution of 1'5, which 
is only slightly worse than that of our submillimeter datasets. 

3. Observational results 

Figure Q] presents maps of integrated intensities of [Ci], 12 CO 
1-0, I3 CO 1-0 and 12 CO 4-3 observed in IC 348. All maps are 
integrated over the velocity range Vlsr = 2 km s to 14 km s _1 . 
The energy of the upper level for the [Ci] 1-0 transition 
is 24 K and the critical density fo r this transition is 10 3 cm~ 3 
for collisions with molecular H2 dSchroder et all 1 199 lb . This 
implies that [Ci] 1-0 is easily excited and the line is eas- 
ily detectable even when emitted by moderate density inter- 
stellar gas exposed to the average interstellar radiation field. 
Previous observations have found the [C 1] to be extended and 
well correlated with the low-7 CO emission (see a review by 
Preibisch et all U996t ISchneider et all 120031 : iMookeriea et al. , 



http://www.iram.fr/IRAMFR/GILDAS/ 
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Fig. 1. Velocity integrated intensities of [C i] emission (color) overlayed with contours of (a) 12 CO 4-3, (b) 12 CO 1-0 and 
(c) 13 CO 1-0 integrated intensities at a common angular resolution of 70". The center of the maps is at a - 03 44 m KH, 5 - 
32°06'(J2000). All tracers are integrated from ~Vlsr 2kms~ I to 14kms _I . Contours range between 20% to 90% with a step of 
10% of the peak intensities that are 73Kkms' for 12 CO 4-3, 72Kkms- 1 for 12 CO 1-0 and 31Kkms-' for 13 CO 1-0. The 
seven filled triangles indicates the positions where we carry out a detailed PDR analysis later. The open pentacle symbol denotes 
the position of HD 28 1159. 



20061) . In IC 348 the [C i] emission peaks to the southwest of 
the mapped region with a value of 24Kkms~ I and we detect 
[C i] emission over a large region with a homogeneous inten- 
sity distribution at a level of 65% of the peak intensity. The 
[C i] emission is extended towards the east and northeast of the 
mapped region. This is consistent with the clumpy UV irradi- 
ated cloud scenario. 

The 12 CO 4-3 emission peaks at almost the same position 
as [Ci] with an intensity of ~ 66Kkms -1 . However, the mor- 
phology of the 12 CO 4-3 emission is different to that of [C i]. 
We find secondary maxima in the north and a steeper decay 
of the emission to the south. The distribution of th e 12 CO 4- 
3 emi ssion agrees very well with that of 12 CO 3-2 (iSun et al 



2006b . We attribute the difference in the intensity distributions 
of [C i] and 12 CO 4-3 to the fact that 12 CO 4-3 emission arises 
from regions of higher temperature and density, while [Ci] 
emission can arise also from embedded PDR surfaces with 
moderate to low density within the molecular clouds. 

The distribution of the 12 CO 1-0 emission shows the same 
global structure as the 12 CO 4-3 emission, but it is considerably 
more extended. This is consistent with the fact that I2 CO 1-0 
traces lower temperature and density. The morphologies of [C i] 
and 13 CO 1-0 appear to be very similar, implying that they 
trace the same material and both lines are effectively column 
density tracers. 

We selected seven positions within IC 348 for a more de- 
tailed study. Six of the seven positions are oriented along a cut 
from the northern edge of the cloud, past HD 281 159 and into 
the clouds. The seventh position is at the southeast edge of the 
cloud. Though emission at the seventh position (southern rim) 
is weak, it is clearly detected. The selected positions vary in 
their physical and chemical conditions. Spectra of [Ci], 12 CO 
4-3, 12 CO 1-0 and I3 CO 1-0 at those seven positions are dis- 
played in Fig. |2] 



The [Ci] emission is centered at a velocity of ~ 8kms _1 . 
The lowest [C i] main beam brightness temperature (~ 2 K) oc- 
curs at the southern rim (9', -10'), while the highest tempera- 
tures of 7.0 K and 6.7 K are observed at the position of the [C i] 
peak (-l',-3') and the western rim (-5', -7'), respectively. The 
line widths of [C i] spectra at the seven positions range between 
2.0kms~ 1 and2.6kms~ 1 . 

The 12 CO 4-3 main beam brightness temperature varies be- 
tween ~ 3.6 K at southern rim and 16.4 K at [C i] peak. The line 
width of the 12 CO 4-3 spectra ranges between 2.0 km s 1 and 
2.7kms _1 . 

Both the 12 CO 1-0 and 13 CO 1-0 spectra consist of two 
velocity components: one is at ~ 7kms _1 and the other lies 
at ~ 8 km s . The two components are most prominently visi- 
ble at the southern rim, map center and the [C i] peak. For the 
8kms _1 component at the [Ci] peak, the 13 CO 1-0 spectrum 
peaks at the dip of the 12 CO 1-0 spectrum, which may indicate 
self-absorption effect for this component. 

Table [Upresents the integrated intensities and line ratios at 
the seven selected positions. Figure[3]shows the variation of the 
integrated intensities and their ratios at these seven positions. 

In Fig. [3^ we see that the intensities of [C i], 13 CO 1-0 and 
12 CO 4-3 lines show similar trends, only the falling off of the 
intensity of 12 CO 4-3 is somewhat less drastic than [C i], 13 CO 
1-0. As opposed to the single peak seen in [C i], 13 CO 1-0 and 
12 CO 4-3, the 12 CO 1-0 intensity profile shows two peaks with 
comparable intensities, one at the northern rim and the other at 
the position of the [C i] peak. 

Three independent line ratios, i.e., [Ci] / 12 CO 4-3, [Ci] / 
13 CO 1-0 and 12 CO 4-3 / 12 CO 1-0, are presented in Table [T] 
and Fig.rjb (in ergs" 1 sr 1 crrT 2 ). The largest [Ci] / 12 CO 4-3 
ratio (1.21) occurs at the southern rim and the second (0.67) and 
third (0.41) largest are at the western rim and [Ci] peak, while 
the ratios at the other four positions are ~ 0.25. When compar- 
ing with other star forming regions like W 3 Main, S 106 and 
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Fig. 2. Spectra of the four different tracers on those six positions along the cut from the northern edge of the cloud into the cloud, 
and on the seventh position at the south of the cloud. [C i] spectra in all the positions have been multiplied by 3. All panels have 
the same x axis from 3.5 km s~ : to 14.5 km s -1 and the same y axis from -0.5 K to 30.5 K. 



Table 1. The observed integrated line intensities (in erg s 
on the integrated intensity are estimated to be ~ 15%. 



sr ' cm 2 ) and line ratios at the seven selected positions. The errors 



Position 


(Aa,AS) 


[Ci] 


12 CO 4-3 
[1(T S ergs" 


12 CO 1-0 

1 -1 -2i 

sr cm -J 


13 CO 1-0 


[ClJl-0 
'-C04-3 


[Cl]l-0 
13 COl-0 


'^04-3 
l-COl-O 



Northern rim 
Cluster 
Ridge 
map Center 
[C i] peak 
Western rim 



( 8', 6') 
(6', 4') 
(4', 2') 

d'.-i') 

(-1V3') 
(-5',-T) 



96.56 

81.99 

75.09 

106.79 

226.61 

172.05 



Southern rim (9' ,-10') 72.71 



336.43 
365.30 
313.20 
413.27 
548.04 
256.15 
60.21 



9.59 
7.14 
5.37 
8.38 
9.31 
6.01 
5.50 



1.74 
1.41 
1.61 
2.10 
3.70 
2.13 
1.35 



0.29 
0.22 
0.24 
0.26 
0.41 
0.67 
1.21 



55.51 
58.01 
46.59 
50.85 
61.23 
80.61 
54.02 



35.07 
51.15 
58.38 
49.33 
58.85 
42.61 
10.95 



Orion A etc. (see Table 3 by iKramer et all |2004 . the [C i] / 
'"CO 4-3 ratio at the other five positions are within the range 
found in those star forming regions except for the western and 
southern rim. The highest ratio at the southern rim is close to 
the ratio foun d at the center of M5 1 and the nucleus position 
of NGC4826 JIsrael&BaasLl2002h . The 12 CO 4-3 / 12 CO 1-0 
ratio peaks at [C i] peak and decreases for the positions further 
away with increasing distance. The [Ci] / 13 CO 1-0 ratios are 
rather constant except for the western rim. 



4. LTE Analysis 

Figure [4] shows a map of the 12 CO 4-3 / 12 CO 1-0 line ratio 
in terms of line integrated temperatures (K km s _1 ) overlaid by 
12 CO 4-3 integrated intensities. The ratio has its minimum of 



~ 0.3 at the edges of the cloud; ratios of ~ 0.9 are found at the 
'"CO 4-3 peaks; maximum ratios between 1.1 and 1.5 occur 
close to HD 281 159. 

For a rough estimate of the column d ensities and for an 
easy comparison to other [Cil studies (e.g. iLittle et all 119941 ; 



Plume et all l2000t Bensch et all 120031) we first apply a Local 



Thermodynamic Equilibrium analysis. Assuming LTE and op- 
tically thick emission, the 12 CO 4-3 / 12 CO 1-0 ratio of 0.3 
corresponds to an excitation temperature T ex of ~ 10 K. The ra- 
tio grows to 0.85 for T ex ~ 50 K and remains almost constant 
at even higher temperatures. Considering the calibration uncer- 
tainty, we set r ex to 50 K in all the positions where the ratio falls 
above 0.85. We note that [C i] and 12 CO often have different ex- 
citation temperatures. However, JV(C)lte only changes by 20% 
when the assumed excitation temperature varies between 20 
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Fig. 3. The top panel (a) presents the normalized integrated in- 
tensities of [Ci], 12 CO 4-3, 12 CO 1-0 and 13 CO 1-0 in the 
seven selected positions. The typical errors of the intensities 
are about 15%, resulting in errors of the line ratios of about 
21%; (b) displays integrated line intensity ratios. The x axes of 
all panels are the relative distance to the northern rim (8', 6') 
position. The southern rim is plotted in a separate box. 
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Fig. 4. The color plot represents the 12 CO 4-3 / 12 CO 1-0 line 
intensity ratio (in units of Kkms -1 ) observed in the IC348 
cloud. The contours show the I2 CO 4-3 line intensity spaced 
at 10% intervals from 20% to 90% relative to the peak inten- 
sity of 73 Kkms -1 . White pixels in the figure denote positions 
where either of the line intensities falls below the 3 <x threshold. 



tio [ 12 CO]/[ 13 CO] of 65 dLangeretallll990l) we translate the 
13 CO column densities into CO column densities. H2 column 
densitie s are derived from th e [CO]/[H2] abundance ratio of 
8 10~ 5 (IFrerking et all Il982l) . The results are summarized in 
Table |2] 

The resulting map of C/CO column density ratios is pre- 
sented in Fig. [5] In most of the cloud, the C/CO ratio falls be- 
low 0. 1 . Higher values (up to 1 .5) occur at the rim of the cloud 
where the 13 CO 1-0 emission is more diffuse. Values as low 
as 0.02 are found at the southwest of the cloud and along the 
ridge to the northeast. As 13 CO 1-0 roughly traces the H2 col- 
umn density, this plot indicates an anti-correlation between the 
C/CO ratio and th e H? column density. T his is consistent with 
the conclusion of iMookeriea et al.l (120061) that the [C 1] line is 
not a tracer of optical extinction, total H2 column density or 
total mass in the beam. 




Aa ['] 

Fig. 5. The C/CO column density ratio overlaid by 13 CO 1-0 
integrated intensities contours. The contour levels are spaced at 
10% intervals from 20% to 90% relative to the peak intensity of 
31 Kkms"'. White pixels in the figure denote positions where 
the 13 CO 1-0 intensities fall below the 3 cr threshold. 



The seven positions, for which we perform a detailed anal- 
ysis within the context of a PDR model, exhibit C/CO column 
density ratios between about 0.05 and 0.5, i.e. they cover a rel- 
atively wide range (see Table |2). The ratios measured within 
the full map cover almost the full range cited so fa r in the liter- 
ature with values as low as 0.03 in parts of G34.3 (Little et al 



1994) and as high as 1 .5 seen in the Polaris Flare (IBensch et al 

2003). 



and 150 K. For the first order estimation, it is therefore rea- 
sonable to use the 12 CO excitation temperature also for [Ci]. 
With the 12 CO excitation temperatures and the assumption of 
optically thin [Ci] and 13 CO 1-0 emission, we then obtain the 
column density of C, JV(C) L te, and 13 CO, JV( 13 CO) L te, from 
the integrated line intensities. Using the relative abundance ra- 



5. PDR Modelling 

In the next step we attempt to obtain a self-consistent model 
for the physical and chemical structure of the cloud matching 
the observed C/CO ratios. We will use a clumpy PDR model 
to self-consistently compute the abundance of the species and 
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Table 2. Results of the LTE analysis. The excitation temperature T ex is listed in Column (2); Column (3) - (5) present C, CO and 
H2 column densities; the C/CO ratio is listed in Column (6). The uncertainties of the column densities are derived by varying the 
integrated intensities by ±15%. 



number 


T 

1 L'X 


AT(C)lte 


/V(CO) LTE 


/ V(H2)lte 


C/COlte 




K 


10 16 [cnr 2 ] 


10 17 [cnr 2 ] 


10 21 [cnr 2 ] 




Northern rim 


18- 


12.08^' 


9 56 2284 


12 50 28 - 55 


13 028 


Cluster 


50 20 


9 61 1117 


16.99^ 4 


21 24 24 - 43 


06° 15 


Ridge 


50 28 


9 04 10 - 46 


19 37 2303 
Ly - Jl 10.43 


74 7 1 28.78 
"""13.04 


05 010 


Map center 


48 50 

^°18 


13 09 1600 


24 39 25 - 03 

■"•-"law 


30 - 49 12'82 


05° 16 


[C 1] peak 


50 29 


26.85 3 ™ 8 


44.47^ 4 


55 59 63 - 93 


06 013 


Western rim 


28- 


19 72 2830 

ly - '^17.61 


16.25f* 


2031 3 ^ 


12 030 


Southern rim 


7* 


37 50 7257 

-"• Ju 21.96 


6.66ft 


8.32^ 7 


56 1 34 



their excitation determining the line strength. The model is con- 
strained by the strength of the UV field, which can be indepen- 
dently estimated, and by the known geometry, provided by the 
complex, filamentary, and clumpy structure of the region. 

5.1. FUV intensity 

5.1 .1 . Estimate from the FIR continuum 

We use HIRES pro cessed 60 and 100 fim IRAS data 
dAumann et all 1 19901) to generate a far-infrared intensity map 
/fir in IC348. It represents the far-infrared intensity between 
42.5 /mi and 122.5 fim as measured by t h e filte r curves of 
the two IRAS data sets (INakagawa et all 119981) . Following 
Kaufman et all (11999b . we assume that the total FUV Aux,;^fir, 



absorbed by the dust grains is re-radiated in the far-infrared. 
In this w ay we estimat e the F UV field from the far-infrared 



field (cf. iRramer et all 1200 5): ^fir/Vo - 4-n /fir, where 



y n - 2.7 x 10~ 3 erg s" 1 cm' 2 dDrainelll978l;lDraine & Bertoldil 



19961) . The derived spatial distribution of the FUV intensity is 



presented in Fig. [6] In the mapped region, the FUV field ;^p IR 
varies between about 1 and 100 Draine units. For the seven se- 
lected positions ^fir is listed in Table [3] 

5.1 .2. Estimate from the stellar radiation 

The primary source of UV radiation in the cloud is HD 28 1 159. 
Assuming that the star is a black body at an effective 
temperature T e g corresponding to its s pectral type, we can 



estimate the FUV flux independently dKramer et al 
2000h 



Beuther et al. 



19961 



A B 5 V star has an effective tempera- 
ture lo g(T P ff) - 4.182 and a total luminos ity L of log(L/L ) = 
2.681 dde Jager & Nieuwenhuiizerlll9871) . The FUV luminos- 
ity Lfuv an d me FUV flux ^fuv are defined as 



tD 



uv 



^FUV 



o„ 



and ^fuv 



^FUV 

4^2' 



(1) 



where <t>uv is the FUV flux between 6 eV and 13.6 eV, O tot ai 
is the total flux, and d is the distance to the UV source. We 
compute the distance d in two ways: for the point close to 
HD 281 159 w e use the opti cal radius of the cluster IC348 of 
about 0.37 pc (IHerbigll 1998b as the minimum distance between 
star and cloud. For all other points we assume that star and 
cloud are located in the same plane so that the distance is di- 



- 20 




a« [■] 



Fig. 6. FUV intensity distribution (contours) in IC 348 esti- 
mated from the IRAS FIR fluxes, overlayed on the [Ci] in- 
tegrated intensities (colors). The FUV intensity ^fir contour 
levels run from 5, 10, 20 to 100 by steps of 20 Draine units, xa- 

rectly given by the observed separation. Finally we normalize 
the FUV field to units of xo- The resulting FUV fields, x^tm, are 
also listed in Table|3]for the seven selected positions. 

We see that both methods provide consistent FUV field val- 
ues. Because of the uncertainty of the distance between star and 
cloud we consider the values ^fir, derived from the HIRES 
data, somewhat more reliable so that we will use them in the 
following. The FUV field x at the seven positions ranges from 
about 1 to 90 Draine units. 

5.2. Clumpy PDR scenarios 

Molecular clouds are known for their highly filamentary and 
clumpy, often self- s imilar structure over a wide range of 



20011: iHever et all 1998; Stutzkietal. , 



scales dDame et al 

1988b . The observed structure further breaks 



structures with every step of increased 



up into sub- 
spatial resolution 



dFalgarone & Phillips! 1 19961: iBensch et all l200lb . For IC 348, 
the self-similar structure was directly measured in terms of 
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the A- variance by ISun et al.l (120061) . In such a cloud, the 
UV field can penetrate deeply into the cloud forming PDRs 
at many internal surfaces. The cloud is basically filled by 
surfaces ( IQssenkopf et all 120071) . This explains the observed 



[Cil and eve n [Cn 



Stutzkietal.. 1988; Howe et al 



Meixner & Tielens, 



emission deep within molecular clouds 



199 It IStacey et all [19931 



19931) . We can simulate such a structure 



by an ensemble of individual clumps, each of them being small 
compared to the overall cloud, leaving enough "e mpty" regions 



to pe rmit a deep penetration of the UV field (IStutzki et al 
1998b. 



In the following we test two different clump ensembles with 
respect to their ability to reproduce the observed line intensi- 
ties. 



5.2.1. Modelling of individual clumps 

All individual clumps in the PDRs are modelled us i ng the 



spherical KOSMA - t PDR model dStorzer et all Il996t 
Rollig et all I2006J I2007I) . The model computes the chemi- 
cal and temperature structure of a clump illuminated by an 
isotropic FUV radiation field and cosmic rays. The emission 
from the model is calculated as a function of the gas density, 
FUV radiation field and mass of the clumps (implicitly speci- 
fying the clump size). The model clumps are assumed to have 
a power-law density profile of n(r) oc nor for 0.2 < r/r c \ < 1 
and n(r) = const, for rjr c \ < 0.2. «o is the density at the clump 
surface called clump density, which is about half of the clump 
averaged density (n): «o ~ 757 («)• 

The PDR models are pre-computed on a regular grid with 
equidistant logarithmic steps. The FUV field x covers the range 
from 10°, 10 05 , ..., 10 5 5 , 10 60 Draine units; The clump densi- 
ties range from 10 20 , 10 25 , ..., 10 55 
masses M c \ cover 10 -30 , 10" 2 ' 5 



10 cm ; and the clump 



10 l .9 10 20 M 



5.2.2. Ensemble of identical clumps (Ensemble 1) 

In a first approach we will assume that the cloud is composed 
of a collection of identical clumps, so that the full ensemble 
is characterized by the parameter set of the individual clumps 
and their number only. To determine the properties of the in- 
dividual clumps we perform a x 1 fit of the three line inten- 
sity ratios given in Table [TJ using a spherical PDR model in a 
three-dimensional parameter space spanned by the FUV field, 
the clump mass, and density. The results for the seven selected 
positions are summarized in Table [3] Figure [7] presents an ex- 
ample of the x 1 fitting for the map center (1', — 1'). The left 
panel shows a cut through the parameter space for a fixed UV 
field, demonstrating the fit of clump mass M c / and clump den- 
sity «o- The right panel presents the corresponding fit for the 
FUV field x an d the clump density «o- We see that all line ra- 
tios are good density tracers, but hardly constrain the clump 
mass as expected. The clump density has to fall between 2.5 
and 5.0 10 5 ctrT 3 . But the reduced^ 2 only changes between 1 
and 2.5 when the clump mass varies from 10~ 20 to 10 10 M Q . 
This general behaviour applies to all positions analyzed. The 
clump densities at the seven positions vary between 4.4 10 4 to 



4.3 10 5 cm 4 . The highest clump density occurs at positions of 
the ridge and the map center. The clump masses are not well 
constrained, the best fitting values fall between 0.1 and 0.4 M . 
The fitted FUV flux ^pdr agrees with ^fir derived from the 
FIR continuum in most positions. The maximum deviations by 
a factor three correspond to the maximum deviation that we 
got comparing the FIR continuum with the stellar radiation es- 
timate. 

The total number of clumps in the beam A^ns is calcu- 
lated by comparing the observed absolute line integrated in- 
tensities with the models, correcting for the beam filling fac- 
tor^! A 70" beam corresponds to a diameter of O.lOpc at the 
distance of IC348. The clump size obtained from the best 
fit models found in the previous section ranges from 0.01 to 
0.02 pc. The beam filling factor ?7bf,i of each clump is cal- 
culated as 77bf,i = Q c i/Qbeam, where Q c i and Qbeam are the 
solid angles of a clump and the beam. Q c i = jtR^/D 2 is the 
solid angle of the individual clump. D is the distance of the 
cloud; assuming a Gaussian beam shape, fibeam i s computed as 
^beam = ^beam/^ m ^' wriere #beam is the FWHM of our beam. 
The intensity correction for a single clump is I! - I m b,Jt]bf,i, 
where 7 m b,i is the beam averaged line intensity for that single 
clump. Since [Ci] is the most optically thin one among our 
four tracers, we use its line intensity to derive the number of 
clumps in the beam iV ens (see Table|4]i by dividing the observed 
intensities /[Ci]obs by the intensities for a single clump /[Cijmod.i 



N a 



'[Cljobs 



I[C rjobs ^bearn 

X 



^[C I]mod,i • 77bf,i I[C I]mod,i ^cl 



(2) 



Then we use N ens to calculate the ratios R of line intensities 
between model ensemble and observations for the other tracers: 



R 



^»ens " MTiodj * ^?bf,i 

^obs 



(3) 



The R values are listed in Table [4] This allows easy judg- 
ment of the quality of the model fit. Considering the calibra- 
tion uncertainty, a perfectly fitting model yields a ratio R of 
1±0.15. All three remaining tracers from Ensemble 1 have a 
good agreement between the modelled and observed absolute 
line intensities. 

5.2.3. Clumps with a mass and size spectrum 
(Ensemble 2) 

Representing the clumpy structure of the cloud by an ensem- 
ble of identical clumps is, of course, an oversimplification. 
A realistic clump size and mass distribution should be taken 
into account. Large-scale CO maps present the clumpy struc- 
ture of the ISM with a mass distribution following a power 
law dN/dM oc M~ a , wh ere N is the number of clumps and 
a has a value around 1.8 (iRramer et al.Ul998bl) . Furthermore, 
the observations also show that there is a strong correlation of 
the density and mass of the clu mps corresponding a mas s-size 
relation r oc M y with y ss 2.3 JHeithausen et all Il998i) . For 



2 By using line integrated intensities, a smaller internal velocity dis- 
tribution in the clumps relative to the total observed line width will be 
automatically corrected by a higher clump number. 
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Fig. 7. Comparison of the observed line intensity ratios [Ci]/ 12 CO 4-3, [Ci]/ 13 CO 1-0 and 12 CO 4-3/ 12 CO 1-0 at (1', -1') with 
the single clump KOSMA - r PDR model calculations. The left figure shows the fit of clump mass and density at a fixed FUV 
field x of 10 L5 ; the right one presents the fit of FUV field and clump density at a fixed clump mass of 10 ' M Q . The central solid 
contours represent the observed intensity ratios, the outer dashed two give those for a 20% uncertainty. The grey-scale images 
indicate the reduced^ 2 of the fit. The position of the minimum reduced^ 2 is marked by a white cross. The horizontal light blue 
line in the right panel shows the FUV Aux^fir from the HIRES data. 

Table 3. Physical parameters at selected positions in Fig Q] obtained by fitting the observed intensity ratios to the single clump 
KOSMA - t PDR models. Columns (2) - (5) list the fitted clump densities, clump mass, clump radius and FUV field; Column (6) 
is the minimum reduced^ 2 of the PDR fit; the FUV field derived from the FIR continuum ^fir and that from the stellar radiation 
^star are in Column (7) and (8); The distance to HD 28 1 159 used for the calculation of ^ sta ris listed in Column (9). 



Position 


n 
[10 5 cnr 3 ] 


M a 

M G 


Rd 
pc 


^PDR 


**• 

■^ mm 


ztFIR 


z^star 


d 
pc 


Northern rim 


1.91 


0.2 


0.01 


15 


1.38 


39 


40 


0.37 


Cluster 


2.75 


0.1 


0.01 


98 


1.05 


84 


40 


0.37 


Ridge 


4.25 


0.4 


0.01 


38 


1.23 


46 


40 


0.37 


Map center 


3.98 


0.1 


0.01 


23 


1.04 


23 


13 


0.64 


[C i] peak 


1.91 


0.1 


0.02 


23 


1.96 


11 


7 


0.90 


Western rim 


0.63 


0.2 


0.02 


20 


2.00 


6 


3 


1.42 


Southern rim 


0.44 


0.4 


0.02 


2 


2.31 


1 


3 


1.29 



an ensemble of randomly positioned clumps with a power law 
mass spectrum, there is a relation among the power law spec- 
tral index (3 of the power spectrum, a and y: /3 - y(3 - a) 
(IStutzki et allll998l) . Using the a and y values above,/? = 2.76, 
which is close t o the powe r law in dex 13 = 2.71 of 13 CO 2-1 in 
IC 348 found bv lSunetalJ J2006I) . 



In this Section we assume these characteristics to be uni- 
versal and model the emission by ensemble averaging the PDR 
single clump result s over such a clump distribution dCubick , 
2005t ICubick et all |2008) to reproduce the observed line in- 
tensities. The remaining parameter space is then given by the 
FUV field input ,Y ens , the total mass within the beam M ens , the 
average clump ensemble density « mea n, and the minimum and 
maximum clump masses, m™ n and »?™ x . 

To constrain the fit to a three-dimensional parameter space 
also for this case, we fix the clump mass limits. The upper mass 
limit has to reflect the total mass of the clump distribution be- 
cause a spectral index a - 1 .8 corresponds to a top-heavy func- 



tion, i.e. a distribution where most of the mass is contained in 
the few most massive clumps. We assume that the total mass 
obtained from the Ensemble 1 model also provides a reasonable 
guess for the total mass in Ensemble 2, consequently setting the 
upper mass limit w™ ax to half of the total mass of Ensemble 1 . 
We set the lower mass limit, m"J m , to 10~ 2 M Q . These small 

' cl ' u 

clumps can onl y exist transiently and evaporate on a time scale 
of ~ 5000 yrs (iRramer et all |2008[). The choice of th e lower 
mass limit is not critical for the fit. lCubick et al.l (120081) studied 
the dependencies of the model output on the clump mass limits 
and found that mainly high - / CO transitions (like 12 CO 8-7) 
are sensitive to the lower clump mass limit. 

The actual fit is performed analogously to the ensemble 
1 case. The results are presented in Table [5] We find a rea- 
sonable agreement between the modelled line intensities from 
Ensemble 2 and the observed ones. The intensity ratios be- 
tween model and observations range between 0.7 and 2.0. The 
fitted FUV field Xeni from Ensemble 2 is consistent with the 
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Table 4. Results of ensemble of identical clumps, Ensemble case 1. The mass and FUV field for each single clump are listed in 
Column (2) and (3); Column (4) - (6) present the line intensity ratios between the ensemble models and observations of 12 CO 
4-3, 12 CO 1-0 and 13 CO 1-0; The C, CO and H2 column densities from the ensembles are presented in Column (7), (8) and (9); 
Column (10) is the C/CO ratios. The total number of the clumps in the ensemble and ensemble mass are listed in Column (11) 
and (12); Column (13) lists the beam filling factor of each clump. The total ensemble mass M ens listed in the table is computed 
as Mens = A^ens x M c \, where M c \ is the mass of a clump. 



Position 


M cI 

M 


Xcnl 


^ 12 C04-3 


^ 12 COl-0 


^ 13 COI-0 


/V(C) cns 
10 16 [cm- 2 ] 


JV(CO) tns 
10 17 [cnr 2 ] 


/V(H 2 U 
10 21 [cnr 2 ] 


C/COens 


AU 


M cns 

M 


^bf.i 


Northern rim 


0.32 


10 


0.63 


0.67 


0.87 


13.97 


18.58 


7.72 


0.075 


4.59 


1.47 


0.17 


Cluster 


0.10 


100 


1.07 


0.99 


1.13 


10.02 


21.56 


8.90 


0.046 


16.75 


1.68 


0.04 


Ridge 


0.32 


32 


0.95 


1.13 


0.95 


8.90 


29.49 


11.29 


0.030 


6.72 


2.15 


0.08 


Map center 


0.10 


32 


0.94 


1.05 


0.97 


13.25 


34.24 


13.05 


0.039 


24.56 


2.46 


0.04 


[C 1] peak 


0.10 


32 


0.94 


1.59 


0.83 


31.54 


29.88 


13.14 


0.106 


24.74 


2.47 


0.08 


Western rim 


0.10 


10 


0.78 


1.43 


0.71 


24.52 


9.33 


6.10 


0.263 


11.51 


1.15 


0.17 


Southern rim 


0.32 


1 


0.84 


0.81 


0.74 


14.62 


6.28 


3.03 


0.233 


5.70 


1.82 


0.17 



Table 5. Results of ensembles with a clump mass distribution and mass-size distribution, Ensemble case 2. Column (2) - (4) 
present the ratios (in erg s -1 sr cm 2 ) of the modelled and observed line intensities for 12 CO 4-3, 12 CO 1-0 and 13 CO 1-0; The 
mean clump ensemble density « me an and the fitted FUV field Xeni are presented in Column (5) and (6). The C, CO and H2 column 
densities from the ensembles are presented in Column (7), (8) and (9); Column (10) is the C/CO ratios. The total masses are 
listed in Column (11); Column (12) is the minimum reduced^ 2 of the fit. 



Position 


«12 C0 4_3 


^ 12 COl-0 


^ 13 COl-0 


[10 5 cnr 3 ] 


*ei,2 


MC)ens 

10 16 [cm- 2 ] 


«(CO) c „ s 

10 I7 [cm- 2 ] 


A'(H 2 ) cns 
10 21 [cm- 2 ] 


C/COe ns 


Mens 

M Q 


x 1 - 

A mm 


Northern rim 


1.00 


1.04 


1.44 


3.2 


3 


14.36 


51.82 


19.28 


0.028 


3.84 


1.43 


Cluster 


1.06 


1.03 


1.13 


3.2 


100 


9.88 


25.58 


9.93 


0.039 


1.97 


0.23 


Ridge 


1.13 


1.26 


0.91 


3.2 


32 


9.05 


23.43 


9.10 


0.039 


1.83 


0.50 


Map center 


1.08 


1.19 


1.08 


3.2 


10 


13.78 


40.71 


15.39 


0.034 


3.07 


0.29 


[C 1] peak 


1.13 


1.73 


0.86 


1.0 


10 


31.24 


40.16 


16.61 


0.078 


2.88 


2.37 


Western rim 


1.83 


2.03 


1.13 


1.0 


10 


30.49 


23.72 


12.61 


0.078 


2.19 


3.37 


Southern rim 


1.06 


0.84 


0.66 


0.3 


1 


13.58 


5.24 


2.62 


0.259 


0.53 


4.91 



field Xem from Ensemble 1 and that derived from the FIR con- 
tinuum in most positions. The maximum deviation occurs at 
the northern rim. 

The reduced^ 2 values are similar for both ensembles at val- 
ues less than 5. In spite of the somewhat more complex model, 
the fit is not better than for an ensemble of identical clumps 
(Ensemble 1). 

To compare the observed absolute line intensities with the 
model prediction, we show two examples of observed and mod- 
elled integrated intensities versus frequency, namely the cool- 
ing curves of 12 CO, 13 CO (upto J - 9-8) and [C 1] for the best- 
fit of both Ensemble 1 and Ensemble 2 models at the northern 
rim and the southern rim (see Fig. [8}. At those two positions the 
modelled and observed line intensities agree very well. Fig. [8] 
also provides a prediction for line intensities of [C 1] 3 P 2 - 3 Pi , 
12 CO 7-6, 13 CO 8-7, etc. Fig.[8]indicates that the tracers stud- 
ied in this paper do not allow to draw any conclusion about 
the clump size spectrum that best fits the observations. Both 
Ensemble 1 and 2 fit the observed data equally well. On the 
other hand, they give very different predictions for the higher J 
CO lines. The line intensities predicted for Ensemble 2 are no- 
ticeably larger than those predicted for Ensemble 1 (see Fig . [8} . 
The 13 CO 8-7 line, tracing the column densities of the warm 
gas, is predicted to be very sensitive to the differences between 
the two models, allowing to discriminate between them. 



Based on the CO and IT column densities from the 
PDR analysis, we find the averaged CO relative abundance 
[CO]/[H2] is ~ 2.4 10" 4 , which i s 3 times larger than th e canon- 



ical ratio [CO]/[H 2 ] of 8 lO -5 dFrerking et all 11982b that we 



used in the LTE analysis. The beam averaged C/CO abundance 
ratio (Table [5]l varies by a factor 10 between 0.03 at the north- 
ern rim, which shows the highest column density of all studied 
positions of 2 10 21 cm 2 , and 0.26 at the southern rim, which 
shows the lowest H2 column density of 3 10 21 cm 2 . The ob- 
served anti-correlation between C/CO and N(Hi) and also the 
absolute values resemble those found in ma ny other Galactic 



molec ular clouds as compiled recently by iMookeriea et al. 



d2006l) . They used the KOSMA-t model to interpret data taken 
in the Cepheus B star forming region. They suggested that the 
emission of high column density peaks is dominated by mas- 
sive clumps exhibiting low C/CO ratios, while positions of low 
column densities are dominated by smaller clumps, exhibiting 
higher C/CO ratios. This scenario is confirmed in IC 348 by a 
more complete analysis using clump mass ensembles following 
the canonical mass and size distributions. However, the scatter 
of the present data is large. Observations of the higher lying 
transitions of [C 1] and CO are needed to further improve the 
analysis. 



10 



Sun et al.: The Photon Dominated Region in the IC 348 molecular cloud 



1E-6 
1E-7 
1E-8 
1E-9 
1E-10 
1E-11 
1E-12 
1E-13 




-Ensemble case 1 
3— Ensemble case 1 
— Ensemble case 1 
Ensemble case 2 
Ensemble case 2: 
Ensemble case 2: 
Observed ,2 CO 1-0 and 4-3 
Observed 13 CO 1-0 
Observed CI 1-0 



CO 
l3 CO 
CI 

"CO 
"CO 
CI 



200 



400 600 

Frequency [GHz] 



1E-7 
1E-8 
1E-9 
1E-10 
1E-11 
1E-12 
1E-13 




— n— Ensemble case 1 

— O— Ensemble case 1 

— — Ensemble case 1 

ffl Ensemble case 2 

® Ensemble case 2 

Ensemble case 2 

■ Observed 12 CO 1-0 and 4-3 

• Observed "CO 1-0 

Observed CI 1-0 



200 



400 600 

Frequency [GHz] 



Fig. 8. Observed (filled symbols) and modelled line intensities of 12 CO, 13 CO, and [C i] from both Ensemble case 1 (drawn lines 
with open symbols) and Ensemble case 2 (dotted lines, open symbols with a cross inside) at the northern rim (left panel) and 
the southern rim (right panel). To better display the data, we artificially shift the [Ci] frequency by +40 GHz for the observed 
intensity. 



6. Summary and Conclusions 

We have presented fully sampled [C i] and 12 CO 4-3 maps of 
the IC 348 molecular cloud, covering a region of 20'x 20'. The 
observed 12 CO 4-3/ 12 CO 1-0 ratios vary between 0.2 and 1.5. 
High 12 CO 4-3/ 12 CO 1-0 ratios occur near the B 5 star, at the 
cloud center and northern edge of the cloud. 

We have estimated the FUV field from the FIR continuum 
intensities obtained from the 60 and 100/iin HIRES/IRAS im- 
ages. The FUV field in the whole observed region ranges be- 
tween 1 and 100 Draine units. We also used HD 281 159, the 
primary source of UV radiation in this region, to estimate the 
FUV field for the seven studied regions and the FUV field 
varies between 3 to 40 Draine units. 

Using a clumpy PDR model based on the spherical 
KOSMA - t PDR model we fitted the observed line intensities 
and intensity ratios at seven selected positions. The observed 
intensity ratios provide a strong constraint on the clump den- 
sities. They fall between 4.410 4 cirT 3 and 4.3 10 5 cm~ 3 . The 
FUV field fitted in the model falls between 2 to 100 Draine 
units, consistent with independent estimates for the FUV field 
derived from the FIR continuum maps by IRAS and from the 
stellar radiation. We found that both an ensemble of identical 
clumps (Ensemble 1) and a model using a clump size and mass 
distribution following the typically observed scaling (Ensemble 
2) produce line intensities which are in good agreement (within 
a factor ~ 2) with the observed intensities. The models confirm 
the anti-correlation between C/CO abundance ratio and hydro- 
gen column den s ity fo und in many regions and explained by 
iMookerjea et al.1 d2006l) . 

The results of this study strengthen the picture of molec- 
ular clouds as highly fragmented, surface dominated objects. 
Although Ensemble 2 is the natural case in this picture, 
Ensemble 1 fits the observed intensities equally well. We pre- 
dicted line intensities for [Ci], CO, and 13 CO transitions for 
both cases, which can be tested by future observations. The 
higher / CO lines are needed to decide on the necessity of the 



clump mass and size distributio ns and to determine the clump 
mass limits dCubick et al.l 120081) . 
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